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Abstract. We present a new investigation of the Cloverleaf (z=2.56) 
based on the combination of archival HST/WFPC2 data, recent IRAM 
CO (7-6) maps and wide field CFHT/FOCAM images. The deep WFPC2 
observation (F814W) shows a significant overdensity of Isiaw ~ 23-25 
galaxies around the Cloverleaf that we interpret as the presence of a dis- 
tant cluster of galaxies along the line of sight. The Cloverleaf is probably 
the result of the lensing effects of a system which includes, in addition to 
a single galaxy, one of the most distant clusters of galaxies ever detected. 
We have modelled the lens accordingly. 

The high-resolution IRAM/CO map provides for the first time the orien- 
tation and the ellipticity of the CO spots induced by the shear component. 
Velocity - positional effects are detected at the 8a level in the CO map. 
A strong limit can then be put on the size, shape and location of the CO 
source around the quasar. The CO source is found to form a disk- or 
ring- like structure orbiting the central engine at ~ lOOkm/s at a radial 
distance of ~ lOOpc, leading to a central mass of ~ 10 9 M possibly m 
the form of a black hole. 



1. Introduction 

The Cloverleaf is the gravitationally lensed image of the quasar H1413+H7 (14 h 
15 m 46 s . 23; 11° 29' 44". J2000.0) at z = 2.558 showing four spots with angular 
separations from 0".77 to 1".36. Since its discovery (Magain et al. 1988), the 
Cloverleaf has been observed spectroscopically and imaged with ground based 
telescopes in various bands from B to I as well as at 3.6 cm with the VLA. An 
early model of the gravitational lens was derived by Kayser et al (1990). The 
main difficulty with this model lies in the fact that it predicts, for a z=1.44 
lens, a mass of ~ 2.5 10 11 h§ M within 0.7" (6 /igQ 1 kpc) radius, which would 
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Figure 1. Image of the Cloverleaf obtained with the IRAM telescope 
at Plateau de Bure. (a) is the total CLEANed image, (c) the CLEANed 
blue-shifted image, (d) the CLEANed red-shifted image and (b) the 
difference between the CLEANed red and CLEANed blue image. The 
CLEANed CO(7-6) maps were obtained with a natural beam of 0.8" x 
0.4" at P.A. 15 deg. They have been restored with a circular 0.5" beam 
for comparison with HST data. Contour spacing is 1.35 mJy/beam for 
la, 2 mJy/beam for lc and Id and 3 mJy/beam for lb, corresponding 
in each case to 2a. 

correspond to a relatively bright normal galaxy: so far, searches in the K band 
of the predicted 'bright' lensing galaxy have been unsuccessful (Lawrence 1996) 
and this fact remains a mystery. 

New data sets are available : the post-COSTAR HST/WFPC2 UV images 
(partly discussed in Turnshek et al. 1997), maps and spectral informations in 
the molecular CO transitions (Barvainis et al 1994, Wilner et al 1995, Barvainis 
et al 1997, Yun et al 1997, Alloin et al 1997) and CCD images of the Cloverleaf 
over a 5' field of view (CFHT/FOCAM archive) allowing a precise astrometry 
of the field obtained through different wavebands. New results obtained from 
these data are discussed in the following sections. Throughout the paper, we 
use Ho= 50 h$o km/s/Mpc, ^o=l an d A = 0. 

2. The IRAM interferometric CO(7-6) map 

Our previous CO(7-6) IRAM interferometer measurements (Alloin et al 1997) 
have been complemented with observations at intermediate baselines. The com- 
bined data lead to the CLEANed integrated map restored with a 0.5" circular 
beam shown in Figure |l]a. In order to search for a velocity gradient, we have first 
derived the spatially integrated line profile, following the procedure discussed in 
Alloin et al (1997). The new CO(7-6) line profile, shown in Figure 2 of Kneib et 
al (1998), exhibits a marked asymmetry with a steep rise on its blue side and a 
slower decrease on its red side. Excluding the central velocity channel (so that 
the split in velocity is symmetric), we have built the blue (-225,-25 km/s) and 
the red (+25, +225 km/s) maps displayed in Figures [j]c and [l]d respectively. The 
difference between the red-shifted and blue-shifted CLEANed maps (Figure |b) 
establishes definitely the presence of a velocity gradient at the 8a level. 
Measurements of the spot characteristics from the CO (7-6) image have been 
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performed (spot flux ratios, sizes and orientations) through a fitting procedure 
in the visibility domain, as explained in Alloin et al (1997). The final parameters 
are provided in Table 1 of Kneib et al (1998) where the spot sizes are intrinsic 
to the image, i.e. deconvolved by the interferometer beam. Although the mea- 
surements have not been corrected for seeing effects (mean seeing estimate of 
the order of 0.2") the spots A, B and C appear to be definitely elongated (see 
Figure [l|b for spot labels) . 

3. HST/WFPC2 data set 

The HST data of the Cloverleaf have been provided by the ESO/ST-ECF Sci- 
ence Archive Facility (Garching). Two data sets obtained with the WFPC2 are 
used here, obtained by Turnshek in 1994 and by Westphal in 1995. Standard 
reduction procedures using IRAF/STSDAS packages have been applied. The 
absolute photometry was obtained using magnitude zero-points given in Holtz- 
mann et al. (1995). Complete information about the final images can be found 
in Kneib et al (1998). 

Regarding the relative astrometry and spot sizes, the PC observations in fil- 
ters F336W, F555W, F702W and F814W provide similar results (see Table 1 in 
Kneib et al 1998). The four spots are stellar-like with a FWHM ~ 1.5 pixel or 
0.068". 

The absolute photometry and the relative intensity ratios (Table 4 in Kneib et 
al 1998) have been computed with the Sextractor software (Bertin & Arnouts 
1996). No PSF fitting was applied as the four different spots are well separated 
on the PC. The important variation of the intensity ratios in U compared to V, 
R and I band could probably be explained by absorption along the line of sight 
by intervening galaxies (HI clouds at redshift ~ 1.7 will act as very efficient ab- 
sorbers in the U band). Turnshek et al (1997) proposed that dust extinction is 
the explanation with a preference for an SMC-like dust extinction at the redshift 
of the quasar. 

Figure ^ shows the deep F814W image, where the presence of numerous 
faint objects over a 40" region in the environment of the Cloverleaf is striking. 
In order to quantify this effect we have computed the object number density, 
<N>, and its dispersion, <tjv, in the magnitude range 1=23 to 25. We have 
estimated this density in different regions across the image: in a 40" diame- 
ter region around the Cloverleaf where the density contrast is clearly visible by 
eye and in various randomly selected areas of similar size. On the whole frame 
we find a mean value <^> w hoie= 45 ± 20 objects/arcmin 2 (1 a). The region 
around the Cloverleaf has <N>ci= 85 ±30 objects/arcmin 2 , with a peak at 130 
objects/arcmin 2 . The bulk of the objects forming the overdensity are red ob- 
jects with R702W — huw ~ 0.9 (from 0.7 to 1.2) and have a small size (< 0.3"). 
While their morphology cannot be obtained from the HST images, their red 
color suggest that we might be dealing with E/S0 galaxies in a high-redshift 
cluster. This cluster could either be linked to the quasar at redshift z~2.56 or 
along the line of sight. Hereafter, we assume this cluster to be at a redshift close 
to that of the narrow absorption systems observed by Turnshek et al (1988), and 
Magain et al (1988): either 1.438, 1.661, 1.87 or 2.07. For the sake of simplicity, 
we shall adopt hereafter a value of 1.7. 
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Figure 2. Image of the HST field around the Cloverleaf. The four 
quasar spots are not resolved but are located at the center of the WF3 
chip (bottom-right). The galaxies detected in the field are overlayed 
with elliptical contours indicating their centroid, orientation and el- 
lipticity. The white contours are iso-number density of galaxies with 
23 < I < 25 (ranging from 30 to 80 galaxies/arcmin 2 ). A significant 
density enhancement is clearly visible around the Cloverleaf, suggesting 
the presence of a distant cluster of galaxies on its line of sight. 

Although the cluster candidate mentioned earlier is acting as an additional lens- 
ing agent, the small angular separation between the four spots implies that part 
of the lensing effect is produced by a galaxy located amid the four spots of the 
quasar. The detection of the lensing galaxy would be of great importance to 
model accurately the mass distribution of the lens. Yet, its presence is still 
elusive (Angonin et al. 1990; Lawrence 1996). We tentatively used the deep- 
est WFPC2/F814W image in order to detect it. The four spots of the quasar 
have been subtracted using the PSF model provided by a nearby star. We con- 
clude that the lensing galaxy must have / > 24. From modelling considerations 
(Bruzual & Chariot 1993), we find that the galaxy should have an absolute 
magnitude larger than M = —21.3 in order not to be detected. Hence, neither 
the visible nor the near infrared data are sufficient to constrain the mass of the 
lensing galaxy. But if we include the likely presence of the distant cluster, it is 
no longer necessary to invoke a large mass for the lensing galaxy since a fraction 
of the convergence and the deflection angle would be contributed by the cluster. 

4. Absolution registration of images in the UV and the mm range 

In order to derive precisely the shear induced by the gravitational lens on an 
extended source in the quasar, it is imperative to register with a high accuracy 
the Cloverleaf image in a waveband corresponding to point-like images of the 
quasar (i.e. R-band corresponding to rest wavelength 1967 A) and in a waveband 
corresponding to an extended source in the quasar (i.e. the CO(7-6) molecular 
emission). The high precision required, better than 0.2", cannot be achieved 
from the HST data because of their limited field of view. Hence, we used a set of 
R and I bands for wide-field exposures of the Cloverleaf field obtained under 0.6" 
seeing conditions at CFHT by Angonin, Vanderriest and Chatzichristou, with 
FOCAM. The absolute astrometry was performed using 6 reference stars across 
the field, from the Cambridge APM database. This allows absolute positioning 
of the brightest spot A within an accuracy of 0.15" (rms). Then, we have used 
the relative astrometry of spots B, C and D, with respect to A, from the HST 
image (relative accuracy ±0.01"). Finally, the IRAM CO(7-6) image is obtained 
within an absolute astrometric accuracy of 0.1" (rms). Therefore, it is possible 
to register both the optical and the millimeter data in an absolute manner 
(Figure [[]), and within an accuracy of 0.15". 
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5. New model of the Cloverleaf gravitational lens 

The modelling of the Cloverleaf lens is based on the minimization algorithm 
described previously in several papers (Kneib et al 1993, 1996). The model in- 
corporates parameters of the lensing potential through a simple analytical repre- 
sentation of the mass distribution, here a truncated elliptical mass distribution 
(Kneib et al 1998). The constraints used for the gravitational lens modelling are 
as follows: 

1. the relative positions of the four quasar spots from the HST and the inten- 
sity ratio taken in the R and I band, the least affected by dust absorption. 

2. the non-detection of a 5th spot, which puts a limit on the size of the lens 
core. 

3. the position of the cluster center as measured from the overdensity of 
galaxies near the Cloverleaf. 

4. both the cluster and the lensing galaxy are assumed to be at z=1.7. 

The relative intensity ratios and the measured shapes of the CO spots are 
used as a test of the model, and provide in addition information on the size and 
geometry of the CO source in the quasar. 

We have computed two types of model: model 1 which includes an individ- 
ual galaxy with a dark halo at z=1.7 and model 2 which considers an individual 
galaxy and a cluster both at z=1.7. These models are not unique but give sim- 
ilar qualitative results. The parameters for the lensing galaxy and the cluster 
component, as well as the full details of the modelling can by found in Kneib et 
al (1998). 

In order to test, at first order, the lens model with the CO(7-6) map, we have 
assumed the CO source to be elliptical with a gaussian profile. We have fitted 
its position, size and ellipticity so that it reproduces the observed CO image. 
The upper limit of its size is provided by the CO elongated spots A and B which 
are close to merging, but still clearly separable. Thus, the model must predict, 
in the image plane, disconnected isocontours of the A and B spots. We found 
a typical size of 460x230 pc (FWHM) for model 1, and 155x110 pc (FWHM) 
for model 2. In the present case, 1" in the source plane translates into 7.68h§Q 
kpc with the chosen cosmology. A summary of the CO modelling is displayed 
in Figure ||. 

The agreement between the expected amplifications and the observed UV flux 
ratios is relatively good. The comparison between the expected amplification 
and the observed ones in CO is not easy. Indeed, differential amplification due 
to the CO source extent and location with respect to the diamond caustic can 
explain the observed differences. The total amplification of the CO emission is 
~ 18 for model 1 and ~ 30 for model 2. Following Barvainis et al (1997), these 
amplification factors translate for model 1 (resp. model 2) to molecular mass 
M(H 2 )= 3. 10 9 M (resp. 2.10 9 M ) and M(HI)= 3. 10 9 M (resp. 2.10 9 
M ). These mass estimates are in good agreement with the dynamical mass 
computed in the next section, provided uncertainties in the inclination and in 
the conversion from M{CO) to M{H2). 
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Figure 3. Results of the lens modelling of the Cloverleaf superim- 
posed on the HST image, (a) is the HST image overlaid with the CO 
observed (-225, +225 km/s); (b) is the CO predicted for model 1, con- 
volved by the interferometer beam, (c) is similar to (b) but for model 
2. (j) is the the CO predicted for model 1, not convolved by the in- 
terferometer beam, (d), (e), (f) is similar to the first row but for the 
blue emission, (g), (h), (i) is similar to the first row but for the red 
emission. The dotted line is the corresponding critical line, (k) gives 
the position of the best fitted CO sources for the blue (dashed), red 
(dotted) and total (solid) emission, for model 1. (1) is similar to (k) 
but for model 2. The central diamond-shape curve (in (b), (e), (h), 
(k) and (1)) is the internal caustic crossed by the lensed CO source at 
redshift z=2.558. 

6. Constraints on the quasar CO source 

In addition to the size of the quasar CO source, 150-460 pc FWHM, as derived 
above, we can obtain some information on its substructure using the velocity 
gradient observed in the CO(7-6) line. The map corresponding to the blue side 
of the line (-225,-25 km/s) shown in Figure |l|c shows that spot C is stronger 
and slightly displaced inwards with respect to its counterpart in the map corre- 
sponding to the red part of the line (+25, +225 km/s), shown in Figure |l]d. A 
hint of this effect can be found in Yun et al (1997) although a close comparison 
of the IRAM results with the OVRO one is hampered by the fact that in the 
OVRO study the spatial resolution is twice lower, the line profile is not shown, 
the velocity interval considered - 145 km/s - is narrower than ours - 200km/s 
- and not positioned precisely with respect to the line center. Using the two 
lens models obtained previously, we optimized the structure of the CO source 
to reproduce separately the red and the blue images (Figure |||). 
For model 1 [resp. model 2], we find a difference of 0.03" (~ 230 pc) [resp. 0.02" 
(~ 150 pc)] between the center of regions emitting the blue and red parts of the 
CO line. We have tested this procedure against uncertainties in the CO/HST 
images registration. A change by 0.05" (half a CO map pixel) has a very minor 
impact on the positions of the source regions emitting the blue and red parts 
of the line. The quasar point-like visible source appears to be almost exactly 
centered between the blue- and red-emitting regions. This is reminiscent of a 
disk- or ring-like structure orbiting the quasar at a radius of ~ 100 pc [resp. 75 
pc] and with a Keplerian velocity of ~ 100 km/s (assuming a 90 deg inclina- 
tion with respect to the plane of the sky), the resulting central mass would be 
~ 10 9 Mq [resp. ~ 7.5 10 8 M Q ]. Elaborating a more sophisticated (realistic) 
model of the CO source and its link with the BAL feature also observed in this 
quasar should be performed in the future. Yet, such a geometry is consistent 
with the molecular torus of the standard AGN model. The resulting differential 
amplification is likely at the origin of the asymmetry observed in the line profile 
(Figure 2 in Kneib et al 1998). 
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7. Summary and outlooks 

This new analysis of the Cloverleaf reveals that this is probably a complex 
lens which includes a lensing galaxy and an additional distant lensing cluster 
of galaxies. The reality of the cluster toward the Cloverleaf has still to be 
confirmed independently. Yet, most of the faint galaxies around the Cloverleaf 
are found in the same magnitude and size ranges, as expected if they indeed 
belonged to a cluster. If the cluster is at a very large distance, the shift of its 
galaxy luminosity function up to higher apparent magnitude would explain why 
the number-density contrast of the cluster with respect to faint field galaxies is 
lowered down to only a 4 a level. 

This interpetation implies that the lensing galaxy may not be very massive 
and consequently may not be very luminous, helping to explain the mystery 
of the lensing galaxy not having been detected so far. The drawback is that, 
despite the constraint that the shapes of the CO spots provide on the orientation 
of the mass density distribution, it mandates sharing the mass between the 
lensing-galaxy and the lensing-cluster which increases the number of possible 
lens configurations. Further data might reveal the position and the shape of the 
light distribution of the lensing galaxy. 

Such information would be useful to improve the mapping of the CO source, as 
emphasized by Alloin et al. (1997). With the present-day data, the CO source 
is found to be a disk- or ring-like structure with typical radius of ~ 100 pc, 
under the lens composite model of a galaxy and a cluster at z=1.7, leading to a 
central ~ 10 9 M object, typical of a massive black-hole. It is amazing to see 
that a disk with such a small intrinsic size can be spatially "resolved" even at 
an angular distance as large as 1.6 K^q Gpc. 

Although this remains to be confirmed independently, the discovery of a distant 
cluster of galaxies on the line of sight to the Cloverleaf is remarkable because it 
reinforces the suspicion that many bright high redshift quasars are magnified by 
cluster-like systems at large distances. This was already reported from analyses 
in the fields of the doubly imaged quasar Q2345+007 (Bonnet et al 1993; Mellier 
et al 1994; van Waerbeke et al 1997), where a cluster candidate is expected to 
be at z~ 0.75 (Pello et al 1996) and of MG2016 where the X-ray emission of 
the intra-cluster gas has been observed and for which the Iron line (from X-ray 
spectroscopy) gives a redshift z~l (Hattori et al 1997). 
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